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ABSTRACT 



A study of aerosol variations in the Central Pacific was conducted utilizing 
NOAA-9 AVHRR data and concurrent shipboard measurements from the NOAA 
R/V OCEANOGRAPHER during the RITS-88 cruise. The transect was conducted 
7 April to 5 May 1988 along longitude 170° W from latitude 50° N to 12° S. Aerosol 
physiochemistry measurements were provided every 1° of latitude. Satellite obser- 
vations of optical depth, Aerosol Particle Size Index (S 12 ), channel 1 (0.63 //m), 
and channel 3 (3.7 /im) low cloud reflectances were analysed. These parameters 
were evaluated during several naturally occurring events, foremost of which were 
the Gobi desert dust storms and the eruption of Kilauea volcano. By comparing 
shipboard and satellite data, satellite retrieval techniques were verified, shipboard 
measurements were expanded to a regional scale, and the relationship between solar 
reflectance and the microphysical properties of clouds was verified. 



TABLE OF CONTENTS 

I. INTRODUCTION 1 

A. MOTIVATION 2 

B. OBJECTIVES 6 

II. RADIATIVE PROCESSES 8 

III. PROCEDURES 14 

A. DATA ACQUISITION/DESCRIPTION 14 

B. DATA PROCESSING/ANALYSIS 16 

IV. RESULTS 19 

A. SYNOPTIC ANALYSIS 19 

B. AEROSOL SOURCES 24 

C. INTERCOMPARISON OF OPTICAL DEPTH AND CONDEN- 
SATION NUCLEI 30 

D. INTERCOMPARISON OF S, 2 RATIO AND NSS-50^" 35 

E. CHANNEL 1 AND CHANNEL 3 LOW CLOUD REFLECTANCE 40 

V. CONCLUSIONS AND RECOMMENDATIONS 46 

LIST OF REFERENCES 49 

INITIAL DISTRIBUTION LIST 51 



LIST OF TABLES 



2.1 AVHRR Channel Wavelengths (Kidwell 1986) 



v 



LIST OF FIGURES 



1.1 Conceptual Diagram of Possible Climate Feedback Loop: the rectan- 

gles are measurable quantities and the ovals are processes linking the 
rectangles 5 

2.1 Channel 1 (.63 /jm) and channel 3 (3.7 /im) reflectance versus mean 

droplet radius, from Mineart (1988) 13 

3.1 RITS-88 Cruise of R/V OCEANOGRAPHER. Wind vectors are in- 
dicated for the Julian Day when data measurements were taken. ... 15 

3.2 Flowchart for Satellite Pixel Processing 17 

4.1 NOGAPS 1200 UTC 10 April 1988 Surface Analysis 20 

4.2 NOGAPS 1200 UTC 10 April 1988 500 mb Analysis 21 

4.3 Number of Clear Observations - Values are within ± 3 days of Ship’s 

Position 23 

4.4 “Kosa” Transport Mechanism, 10-14 April 1988, from Murayama (1988) 25 

4.5 “Kosa” Transport Mechanism, 16-20 April 1988, from Murayama (1988) 25 

4.6 Optical Depth Composite. Data taken from NOAA-9 AVHRR during 

April 1988 27 

4.7 S 12 Composite. Data compiled from NOAA-9 AVHRR during April 

1988 29 

4.8 Optical Depth - Measured by NOAA-9 AVHRR within ± 3 days of 

Ship’s Position 33 

4.9 Condensation Nuclei - Measured by the R/V OCEANOGRAPHER . 34 

4.10 S \2 Ratio - Measured by NOAA-9 AVHRR within ± 3 days of Ship’s 

Position 38 



vi 



4.11 NSS-SC^ - Data - Measurements made by R/V OCEANOGRAPHER 

during April 1988 39 

4.12 Channel 1 Low Cloud Reflectance - Measured by NOAA-9 AVHRR . 42 

4.13 Composite Optical Depth - Measured by NOAA-9 AVHRR 43 

4.14 Composite Su - Measured by NOAA-9 AVHRR 44 

4.15 Channel 3 Low Cloud Reflectance - Measured by NOAA-9 AVHRR . 45 



vii 



ACKNOWLEDGMENTS 



I would like to thank Professor Philip Durkee for his assistance, guidance, 
and support throughout the course of this research. A special thanks to Professor 
Robert Haney for his careful review of the manuscript. Finally, the staff of the 
Naval Postgraduate School’s Interactive Digital Environmental Analysis (IDEA) 
Laboratory is commended for their support, especially Mr. Craig Motell. Without 
his programming skills, this project could not have been completed. 



viii 



I. INTRODUCTION 



The advent of high resolution satellite sensors have allowed for global climate 
processes to be studied in much greater detail in recent years. Atmospheric elements 
such as suspended aerosols and various types of clouds interact with the electromag- 
netic energy within each wavelength band creating unique radiation patterns which 
the satellite can detect. The scattering and absorption characteristics of these ele- 
ments are responsible for distinctive radiance signatures. Consequently, information 
about identifiable atmospheric features may be derived from careful examination of 
these radiances. With the development of calibrated, multi-channel sensors in con- 
junction with the great increase in computer power, it is now possible to utilize 
objective techniques to study the earth’s environment with satellites. 

Direct observations by ship or aircraft have provided almost the entire body of 
knowledge concerning aerosol particle distribution over the ocean. Although these 
data are accurate and useful, it’s limited area coverage has been a great problem. 
Additionally, this data collection process does not allow for long term continuous 
coverage of a specific area or give the broad coverage necessary for global climate 
studies. 

Satellite remote sensing is currently the only method of continuously moni- 
toring aerosol particles on a global scale. Frost (1988) has developed an algorithm 
which allows the user to detect aerosol particle characteristics from satellite mea- 
surements. Although not able to detect the type of aerosol, measurements can infer 
the relative size (large or small) of marine aerosol particles in a given volume of 
air. When correlated with other in-situ data, a good measurement of the spatial 
and temporal distribution of atmospheric particles occur, and allow for verification 
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of satellite retrieval techniques. This research will perform such a comparison uti- 
lizing aerosol physiochemistry measurements from the Radiatively Important Trace 
Species (RITS)-88 cruise (Clarke 1989) versus satellite-band NOAA-9 AVHRR mea- 
surements for the same time period. 

A. MOTIVATION 

There are two major motivations for this type of data analysis and the exam- 
ination of the multi-spectral radiative characteristics of both cloud-free and cloud 
conditions in the atmosphere. They are separated into a military application and a 
climate application. 

From the military perspective, infrared wavelength image and ranging systems, 
laser-guided weapons systems, and highly advanced laser communication systems 
have all demonstrated extreme sensitivity to environmental factors. The molecular 
and aerosol components of the atmosphere may absorb emitted energy, or scatter 
a large percentage of it, thereby greatly reducing a military electro-optic systems 
effectiveness. Milton (1977) worked with infrared wavelength imaging and ranging 
systems and demonstrated that they are strongly affected by variable distributions 
of atmospheric aerosols. Shipboard laser-guided weapons with their low level tra- 
jectories are especially susceptible to the effects of suspended aerosols. Bloember- 
gen et al., (1987) demonstrated that the optimum wavelength region for any laser 
system is 0. 4-1.0 urn. This wavelength band avoids severe attenuation from atmo- 
spheric molecular components at short wavelengths and w r ater vapor absorption at 
longer wavelengths. However, what must be accounted for in this band is atmo- 
spheric aerosols, as they are strong contributors to energy losses. Cordray et ah, 
(1977) conducted operational tests propagating a deuterium fluoride chemical laser 



2 



through the lower atmosphere at 3.2 //m. The results conclude that aerosol den- 
sity, size distribution, and chemical and physical composition constitute most of the 
parameters affecting laser sensitivity. The implications of these results are that real- 
time microscale measurements of the atmospheric aerosol content is critical for the 
theater commander to assess and predict the performance of electro-optical systems. 
Knowledge of the characteristics and locations of aerosol features is vital to deter- 
mining the sensor settings, weapon paths, and launch/no-launch criteria. Clearly 
the availability of these data would greatly enhance a mission’s success percentage. 

From the climatic perspective, the scattering and absorption of solar radiation 
by aerosol particles and clouds influences the earth’s radiation budget and thus can 
have a significant impact on climate variations. Various hypotheses have been pro- 
posed describing both a global warming and cooling trend. The “greenhouse effect”, 
absorption by the atmosphere of infrared energy emitted from the earth, has been the 
most publicized theory for a global warming trend. Additionally, deforestation and 
volcanic eruptions have been proposed as planetary sources for warming/cooling. 
Twomey et al., (1984) suggests that increased atmospheric pollution may lead to a 
general cooling trend due to increases in small aerosol particles in the atmosphere. 
Increases in aerosol concentration lead to increases in the number of cloud concentra- 
tion nuclei (CCN). Assuming a constant liquid water content (LWC), the increase in 
CCN increases the number of cloud droplets and decreases the droplet radius. This 
results in increased reflection of incoming short-wave radiation and thus a global 
cooling trend. At the same time, high cirrus clouds have the ability to increase 
absorption of outgoing terrestrial long-wave radiation, thereby warming the earth. 
Thus, an important climatic dilemma exists; clouds both reflect incoming short-wave 
radiation (cooling effect) and absorb outgoing long-wave radiation (warming effect). 
Atmospheric temperature fluctuations indicate that many processes are operating 
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to produce the observed mean annual temperature. Coakley et al., (1983) showed 
that the effect of aerosols on the earth’s radiation budget, through their influence on 
clouds, may be several times that of the direct interaction of the aerosol with solar 
radiation. Due to the extreme complexity and spatial variability of clouds, most 
climate models tend to neglect or over-simplify intricate environmental processes 
which tend to control climate. 

One such important process is described by Charlson et ah, (1987) and is 
represented by Fig. 1.1. He proposes that the majority of CCN formed in remote, 
marine regions are due to aqueous dimethylsulfide (DMS) emissions through a series 
of complex biological, chemical, and physical processes. The rectangles in Fig. 1.1 
represent measurable quantities, while the ovals represent processes that may or 
may not be fully understood. DMS is excreted by marine phytoplankton found 
throughout the planet's aquatic environment. This DMS is then ventilated to the 
atmosphere via an air-sea transport mechanism where it undergoes an oxidation 
process. The result of this chemical reaction is the formation of two products: a 
methane sulfonate (MSA) and a sulfate. The sulfate is thought to be the main 
product and undergoes a transformation to non sea-salt sulfate ( NSS — S0 \~ ) 
particles. 

The abundance of these aerosols in the marine layer support Charlson’s hy- 
pothesis that these are the principal CCN particles. Observations by Pruppacher 
and Klett (1978) demonstrated that observable CCN concentrations vary from 30 to 
200 cm -3 whereas the sea-salt particle concentrations are usually not greater than 
1 cm -3 . Therefore, sea-salt particles cannot be the main CCN source. Two possi- 
bilities now exist for the NSS — SOl~ . First, they may add mass to the existing 
CCN which would tend to increase the particle radii. Conversely, these particles 



4 



Number concentration 


^^''Scatterlng of Solaris. 


Cloud 


of cloud droplets 


radiation by droplets ) 


albedo 


(Fixed LWP) 










Figure 1.1: Conceptual Diagram of Possible Climate Feedback Loop: 
the rectangles are measurable quantities and the ovals are processes 

linking the rectangles. 
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could create new and hence smaller CCN. This resultant change in size distribution 
would then affect the reflectance characteristics of the cloud. Twomey (1977) showed 
that this decrease in droplet radius leads to an increase in the total surface area and 
thus an increase in the cloud albedo. This implies that oceanic biological variability 
may play a major role in marine cloud albedo changes and thus in planetary climatic 
trends. 

Identification of variable aerosol contributions to the radiative properties of 
cloud fields, such as continental pollution sources, volcanic sources, and DMS pro- 
duction is a major motivation for conducting this research. 

B. OBJECTIVES 

Satellite measurements of aerosol particles can provide a better understand- 
ing of the complex interrelationship between the atmosphere and the ocean. The 
primary objective of this thesis is to provide a verification of satellite retrieval tech- 
niques with regards to aerosol particle characteristics. Pfeil (19S6) with modifica- 
tions from Frost (19SS) utilized the satellite measured red and near infrared radiance 
to infer the size relationship of aerosol particles. This is accomplished by forming a 
ratio of the albedos, as determined from these two wavelengths. 

The second objective is to utilize the ship-based observations recorded during 
the RITS-8S cruise and expand the data to regional scales utilizing the NOAA-9 
AVHRR satellite in the visible, near infrared and infrared portions of the electro- 
magnetic spectrum. The two primary areas of interest are at 20° N and 40 0 N, where 
large concentrations of aerosol particles exist. The identification of the sources for 
these aerosols will be shown using satellite retrieval methodology. 

Finally, the relationship between aerosol characteristics and cloud reflectance 
properties will be verified. Mineart (1988) demonstrated the ability to deduce cloud 
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microphysical parameters (primarily mean cloud droplet size distribution) of marine 
stratocumulus clouds utilizing AVHRR satellite data. A multispectral image pro- 
cessing technique is used which combines AVHRR channels 1 (0.63 pm), 3 (3.7 pm), 
and 4 (10.9 pm) to obtain reflectance values. These reflectance values, obtained 
from the same region as the RITS-88 cruise data, will allow for a determination 
of the relationship between solar reflectance and the microphysical properties of 
stratocumulus clouds. 
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II. RADIATIVE PROCESSES 



Radiant energy propagating through the earth’s atmosphere is attenuated by 
interactions with gaseous molecules, suspended aerosols, and the cloud’s and earth’s 
surface. This attenuation, which includes absorption and scattering, is highly de- 
pendent on such factors as electromagnetic wavelength and the composition of the 
atmosphere. The AVHRR sensor located on the NOAA-9 satellite measures wave- 
lengths in five channels, each of which is located in an atmospheric absorption 
“window”, where the incoming electromagnetic energy is primarily affected only by 
scattering. Table 2.1 gives a complete listing of the channels and respective wave- 
lengths. Channel 5 thermal emittance measurements will not be utilized in this 
thesis. 

Backscattered and emitted radiance measured by the AVHRR can provide sig- 
nificant information about atmospheric and surface features. Because each channel 
detects radiant energy in a separate wavelength band, different characteristics of 
these features are revealed based on scattering properties of aerosols. Suspended 
aerosol quantities and size distributions along with several cloud characteristics, 
primarily cloud thickness and liquid water content (LWC), are revealed by the mea- 
surements of channel 1 and channel 2 (Frost 1988). Channel 3 radiance measure- 
ments provide different cloud properties, such as the number of droplets and their 
size distribution, which may be inferred from cloud concentration nuclei (CCN) 
characteristics. Channel 4 radiance measurements provide an indication of surface 
and cloud-top temperatures. The analysis of these four channels provide powerful 
tools for studying aerosol and cloud parameters. 
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TABLE 2.1: AVHRR Channel Wavelengths (Kidwell 1986) 



CHANNEL 


RADIANCE 


WAVELENGTH BAND 


1 


Visible 


0.58 


- 0.68 pm 


2 


Red visible/ Near infrared 


0.70 


- 1.10 pm 


3 


Near infrared 


3.50 


- 3.90 pm 


4 


Thermal infrared 


10.30 


- 11.50 pm 


5 


Thermal infrared 


11.30 


- 12.40 pm 



Optical depth is the first aerosol parameter to be discussed. Because satel- 
lites detect upward scattered radiation from aerosol particles, as upwelled radiance, 
Durkee et al., (1986) showed that there is a direct correlation between the upwelled 
radiance measured at the satellite and the aerosol particle optical depth. When the 
optical depth is linked to the size and distribution of aerosol particles, it is possi- 
ble to relate the satellite detected radiance to the size distributions of the marine 
aerosols. Optical depth is defined as • ' 

6= f" f3 ext dz . (2.1) 

Jo 

This is the extinction coefficient of the entire atmosphere integrated from the surface 
to the satellite. The extinction coefficient is given by 

roo 

ft ext ~ / Trr 2 Q(m,r) ext n(r)dr , (2.2) 

where nr 2 is the cross-sectional area of the particle and n{r)dr is the size distribution 
of the particle at the cross-sectional area. Q ext is an extinction efficiency factor that 
varies with the radius (r) and composition through the index of refraction (m). 
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Equations (2.1) and (2.2) show the relationship between optical depth and the 
aerosol particle. f3 ext is a function of aerosol particle size and distribution. Because 
optical depth is a function of fl ext , it follows that changes in aerosol particle size or 
distribution will be reflected in optical depth and satellite-detected radiance. 

When aerosol particles change, there also must be changes in the radiative 
properties of the atmosphere. Durkee (1986) clearly shows that any change in the 
optical depth results in a change in the diffuse upward radiance at the top of the 
atmosphere. Additionally, optical depth is a function of the single scattering albedo, 
satellite zenith angle, solar radiation, and scattering phase function. Based upon 
the previous assumptions, satellite zenith angle and phase function are the major 
variables in determining optical depth. Durkee (1984) showed that satellite zenith 
angle is the most critical factor affecting the radiance-optical depth relationship for 
small scale analysis. Frost (1988) proved that a variable phase function, such as 
the Heney-Greenstein scattering phase function, should be utilized on large scale 
optical depth analysis where large changes in aerosol particle distributions are ex- 
pected. Because optical depth is calculated with a phase function that varies with 
the dominant aerosol particle, a more accurate representation is obtained. 

The second aerosol parameter to be discussed is the relationship between ra- 
diance and the aerosol size distribution. Durkee (1984) showed that the aerosol 
particle dependent terms in the radiative transfer equations (RTE) are larger at red 
wavelengths (channel 1) than at near infrared wavelengths (channel 2). Because dif- 
fuse radiance is proportional to aerosol particle dependent terms, spectral variations 
in radiance result from variations in aerosol particle dependent terms. This ratio, 
labeled the aerosol particle size index (£ 12 ) is quantified, as 

C )red ) re d 

12 “ (L A ) nir ~ (Mn, r ’ 
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where 6a is the optical depth due to aerosol particles. The ratio of channel 1 to 
channel 2 reflectance enhances the variation between aerosol size distributions in 
the cloud-free atmosphere. More small particles will display a significantly greater 
difference between channel 1 and channel 2 than a size distribution weighted toward 
more large particles. Although satellite derived measurements cannot provide the 
actual particle size distribution curve, the greater the Su ratio, the greater the 
influence of small particles relative to large particles in the atmosphere and vice- 
versa. 

The final aerosol parameter to be discussed is the change in the cloud re- 
flectance properties based on CCN/aerosol variations. To evaluate these changes, 
it is necessary to evaluate both the AVHRR channel 1 and channel 3 reflectance. 
For channel 1, since there is no absorption at this wavelength, incoming photons 
are either transmitted through the cloud or scattered by cloud droplets and re- 
flected. If the cloud thickness is increased, there is less likelihood that energy will 
be transmitted through the cloud, thus increasing backscattered radiation. If the 
LWC is increased, backscattered energy will also increase due to more interactions 
between photons and cloud droplets. Therefore, channel 1 reflectance is a function 
of size distribution, cloud thickness, and LWC. This makes it very difficult, to assess 
variations in reflectance values without making limiting assumptions. 

Conversely, for channel 3 reflectance, the dependence on LWC and cloud thick- 
ness is not nearly as great as the dependence on droplet radius. The sensitivity of 
reflectance to droplet radius in channel 3 is greatest, due to the moderate absorption 
at 3.7 pm (Coakley and Davies 1986). Therefore, there is virtually no transmission 
through the cloud at this wavelength. This means that LWC and cloud thick- 
ness variations do not appreciably affect the cloud reflectance. This allows one to 
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conclude that, at channel 3 wavelengths, the reflectance would increase based on 
injection of smaller CCN which leads to smaller and more numerous cloud droplets. 

Mineart (1988) further illustrated this point by comparing strato-cumulus 
cloud reflectance to mean droplet radius at both channel 1 and channel 3 wave- 
lengths. Fig. 2.1 shows the results of these measurements. There is no correlation 
of channel 1 with droplet size because of the competing effects of size distribution 
with cloud thickness and LWC. For channel 3 reflectance values, as expected, a 
strong linear correlation occurs. As the mean droplet radius increased, the channel 
3 reflectance values decreased, indicating that changes in cloud reflectance based on 
size distribution are observable in channel 3. 
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Figure 2.1: Channel 1 (.63 //m) and channel 3 (3.7 //m) reflectance versus 
mean droplet radius, from Mineart (1988) 
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III. PROCEDURES 



A. DATA ACQUISITION/DESCRIPTION 

The satellite data used in this thesis were obtained from the National Environ- 
mental Satellite Data and Information Service (NESDIS) digital archives. The data 
were collected by NOAA-9, a sun-synchronous near polar orbiting satellite with a 
nominal altitude of 833 km. Upwelled radiance was measured by the Advanced Very 
High Resolution Radiometer (AVHRR) sensor on the satellite. Maxjmum resolution 
of data at the satellite sub-point is 1.1 km x 1.1 km. Global Area Coverage (GAC) 
data were examined during the period 7 April - 5 May 1988. The use of GAC data, 
which minimizes archiving requirements, reduces the resolution to approximately 4.0 
km x 4.0 km. The primary area of concern for this thesis is a remote marine region 
bounded by 50° North and 20° South latitude and 175° and 165° West longitude. 

The Shipboard measurements were obtained from the R/V OCEANOGRA- 
PHER, a National Oceanic Atmospheric Administration (NOAA) funded research 
vessel, during the Radiatively Important Trace Species (RITS)-88 cruise conducted 
7 April 1988 to 5 May 1988. Sampling began at 50° N latitude and terminated 
at 12° S latitude along longitude 170° W. A complete ship’s track is displayed in 
Fig. 3.1 along with wind observations. Data measurements were obtained for every 
1° of latitude and consisted of Dimethyl Sulfide (DMS) emissions in the sea and 
air, NSS — SOl~ aerosol concentration, aerosol size distribution, and total particle 
concentration (Bates and Charleson 1989). 
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CRUISE TRACK RITS 1988 




Figure 3.1: RITS-88 Cruise of R/V OCEANOGRAPHER. Wind vectors 
are indicated for the Julian Day when data measurements were taken. 
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B. DATA PROCESSING/ANALYSIS 

Processing was performed at the Naval Postgraduate School’s Interactive Dig- 
ital Environmental Analysis (IDEA) Laboratory located in Monterey, California. 
The AVHRR satellite data were processed in order to derive the desired reflective 
and thermal parameters and to produce the associated images. Each of the six 
NOAA-9 GAC tapes (72 individual satellite passes) processed required approxi- 
mately 19 hours of CPU time per tape to obtain the necessary data to cover the 
entire area of the RITS-88 cruise. 

The AVHRR data were analyzed based on a scheme developed by Pfeil (1986) 
and modified by Frost (1988). Fig. 3.2 is the flowchart by which the data were 
processed. The applicable data sets produced from this algorithm are as follows: 

1. Channel 1/Channel 2 radiance ratio for cloud-free atmosphere. 

2. Channel 1 reflectance for low clouds. 

3. Channel 3 reflectance for low clouds. 

4. Optical depth using a variable Heney-Greenstein phase function for cloud-free 

atmosphere. 

5. Number of partly cloudy pixels. 

6. Number of high cloud pixels. 

Results for each data set were stored and averaged in 1° x 1° boxes. The 
number of pixels in each box varied depending upon the series of pixel tests in 
accordance with Fig. 3.2. 
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Figure 3.2: Flowchart for Satellite Pixel Processing 
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The first test in the algorithm averaging process is to verify that the input 
AVHRR pixel is valid. This is done by determining if the channel 2 or channel 4 
radiance values are greater than zero after Rayleigh radiances have been removed. 
If less than zero, the input pixel is not processed. The next test is to check for 
contamination due to specular reflection. Because sunglint is a function of satellite 
viewing geometry, it is removed by calculating an elliptical shaped viewing area. If 
sunglint contamination is present, the pixel is not processed. The third test in the 
algorithm is to determine if the pixel is contaminated by high clouds. Should the 
channel 4 brightness temperature be less than 273 K, then the pixel is considered to 
be high cloud contaminated and summed in the high cloud data set. The next test 
performed is a check for low cloud contamination. If the channel 2 cloud albedo is 
greater than 40%, it is considered a low cloud pixel and summed in the channel 1 
and channel 3 reflectance data sets. The fifth test in the series is a preliminary check 
to identify clear pixels. If the channel 1 /channel 2 (5n) albedo ratio is greater than 
1.5, then the pixel is tentatively classified as clear, and is subjected to the smooth 
test for final determination. If not, it is summed in the partly cloudy pixel data set. 

The smooth test consists of finding the maximum and minimum channel 2 
albedo raw counts for the four pixels surrounding the pixel in question. If the 
difference between the maximum and minimum pixels is less than or equal to 5, the 
pixel passes the smooth test, and is summed in the optical depth and (Sn) data sets. 
Otherwise, the pixel is flagged as a partly cloudy pixel, and the process continues 
until all pixels on the data tape have been accounted for. From these data sets, the 
final 1° x 1° images are produced for comparison against the surface shipboard data 
collected by the R/V OCEANOGRAPHER. 



18 



IV. RESULTS 



This section describes the results obtained from the satellite measurements 
and the shipboard observations. First, a brief introduction to the synoptic weather 
situation, major environmental factors, and color-enhanced composite images for 
optical depth and Aerosol Particle Size Index (APSI) for the entire region will be 
discussed. Second, comparisons between shipboard data and satellite data will be 
presented. The satellite data utilized is within a ± 5° longitude band and ± 3 days 
of ship’s position during data collection. The comparison discussions are: 

— Satellite derived optical depth measurements compared to the shipboard 
measurements of condensation nuclei (CN) 

— The Aerosol Particle Size Index (channel 1 /channel 2 ratio) compared to 
the shipboard measurements of atmospheric NSS-^O^ - concentrations 

- Channel 1 and Channel 3 low cloud reflectance values for the entire ob- 
servation period 

A. SYNOPTIC ANALYSIS 

Understanding the synoptic weather pattern is critical to determining the 
source of aerosol particles in this region. Identifying this source helps to explain 
the results obtained from both the satellite and shipboard data. Fig. 4.1 and 4.2 
display a representative surface and 500 mb synoptic flow pattern for the April 1988 
period. 
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Figure 4.1: NOGAPS 1200 UTC 10 April 1988 Surface Analysis 
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Figure 4.2: NOGAPS 1200 UTC 10 April 1988 500 mb Analysis 
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One of the important features affecting the aerosol concentrations in this re- 
gion is the strong westerly flow pattern off the Asian continent, which provides a 
mechanism for the transport of continental and man-made aerosols. Another sig- 
nificant feature is the strong low pressure system centered around 50° N and 170° 
W. This storm system greatly inhibited the shipboard measurements in the 40° to 
50° N region and limited the amount of data available for verification against satel- 
lite measurements. Fig. 4.3 shows the number of clear observations obtained by 
the satellite within a ± 3 day period of the ship’s position. The choice of a 3-day 
band provides a complete data set for comparison with shipboard observations. If 
a shorter time period was utilized, data gaps would exist in the region of 38-44° N 
due to the strong cyclone. Choice of a longer period would result in the satellite 
data less accurately reflecting the in-situ data measured at the surface. 
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Figure 4.3: Number of Clear Observations - Values are within ± 3 days 

of Ship’s Position 
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The last significant synoptic feature to identify is in the region south of 20° N 
where NE trade winds dominate during this period. In the region of the Hawaiian 
Islands, this is especially important as it provides the primary transport mechanism 
for the volcanic emissions from Kilauea (19° N, 157° W). 

B. AEROSOL SOURCES 

In conjunction with the synoptic situation, another key to the study of aerosol 
variations are the environmental influences. Two major influences are present in 
this region during the April 1988 period. These are desert dust storms and volcanic 
emissions. 

“Kosa” is the phenomenon in which dust clouds originating from storms in the 
Gobi desert transport large quantities of dust through the atmosphere. Murayama 
(1987) showed that April is the highest frequency of these dust storms. During the 
period of this study, two significant “Kosa” occurred, which could, when coupled 
with the synoptic situation, affect aerosol measurements in the central North Pacific 
Region. Shaw (1980) has shown that with a nominal zonal translation of 1000 
km/day, dust from storms in the Gobi desert could easily reach the central North 
Pacific. Therefore, these dust storms should have a major impact on aerosol particle 
concentration in the region studied. 

Fig. 4.4 and 4.5 show the four-day transport of particles from these storms 
measured by the geostationary meteorological satellite (GMS). In Fig. 4.4, the dust 
storm originates on 10 April 1988 at 40° N, 140° E. In Fig. 4.5, a similar dust storm 
occurs originating on 16 April 1988 at 42° N, 105° E and results in atmospheric dust 
transport to a region centered at 38° N, 132° E. 
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Figure 4.4: “Kosa” Transport Mechanism, 10-14 April 1988, from 

Murayama (1988) 




Figure 4.5: “Kosa” Transport Mechanism, 16-20 April 1988, from 

Murayama (1988) 



25 



The second major aerosol producing event is the eruption of the Kilauea vol- 
cano on the island of Hawaii (19° N, 157° W). Debris from this eruption was carried 
by the NE trade winds across the track of R/V OCEANOGRAPHER and should 
have a significant effect on the aerosol measurements and allow for verification of 
satellite retrieval techniques and aerosol transport mechanisms. 

Fig. 4.6 is the aerosol optical depth composite for the region from latitude 50° 
N to 12° S and from longitude 140° W to 160° E for the period 7 April to 5 May 
1988. 
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Figure 4.6: Optical Depth Composite. Data taken from NOAA-9 

AVHRR during April 1988 
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In this figure, there are consistently high values for optical depth in the 0.1 to 0.2 
range along the 35° to 45° N latitude band. This is indicative of the transport of 
continental sources by the predominant synoptic flow pattern in the western North 
Pacific. The second key feature is near the Hawaiian Islands depicted by an * in the 
figure. In the 18° to 20° N latitude band, the optical depth increases and spreads 
out as the volcano emissions are transported southwest by the synoptic flow pattern. 
This indicates that either the particles are increasing in size as they move away from 
the source or variations in the rate of emission may be occurring. Fig. 4.7 is the 
Si 2 composite for the same region and time period as the optical depth composite. 
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Figure 4.7: S n Composite. Data compiled from NOAA-9 AVHRR 

during April 1988 
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In the 35-45° N band high S 12 values are collocated with the high optical 
depth values. This is indicative of a very large number of aerosol particles in the 
atmosphere. The high 5 12 values in this same region mean that there are scattering 
agents present. The high 5 12 values in this same region show that the distribution 
of the aerosol particles are weighted towards smaller particles. Therefore, it is 
reasonable to conclude that large amounts of aerosols, continental in origin, are 
blowing off the land and out over water into this region. In the 18-20° N latitude 
band near the Hawaiian Islands, Sn values are decreasing as the volcanic emissions 
spread westward. This indicates that the ratio of small particles to large particles 
is decreasing. When combined with the results of the optical depth increase for 
the same area, it appears that the aerosol particles are increasing in size as they 
transport westward through the atmosphere. 

C. INTERCOMPARISON OF OPTICAL DEPTH AND CONDENSA- 
TION NUCLEI 

Optical depth measured from a satellite gives an understanding of the trans- 
mittance capabilities of the entire atmospheric layer, whereas shipboard measure- 
ments are confined to the surface boundary layer. Nevertheless, a strong indicator 
of the validity of satellite retrieval techniques should be obtained based upon a pos- 
itive correlation between shipboard and satellite measurements. Previously, it was 
shown that changes in aerosol particle distribution detected by the satellite are a 
function of the optical depth. In-situ measurements of condensation nuclei are a 
function of the surface boundary layer distribution of aerosol particles. Therefore, 
if satellite retrieval methods are correct, a positive correlation should exist between 
optical depth and condensation nuclei measurements to the extent that the surface 
reflects the tropospheric changes. 
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Fig. 4.8 shows the optical depths obtained from the satellite measurements 
utilizing a ± 3 day variation from the ship’s position. Optical depth maximums 
occur at 45° N to 35° N, 20° N, 10° N, and 2° N. Additionally, optical depth values 
are minimal in the region south of 2° N, with satellite data not available south of 6° 
S. Fig. 4.9 shows the condensation nuclei count for the same period measured on the 
R/V OCEANOGRAPHER. Similar peaks are observed as with the optical depth 
values. However, there is some variation in their location. There are several reasons 
for these variations in peak location. The primary reason is that satellite measure- 
ments are integrated over the entire atmosphere, whereas condensation nuclei are 
counted at the surface. Therefore, aerosol particles transported into the region from 
distant continental sources should cause a higher optical depth. The second reason 
is that the condensation nuclei count is the sum of all particles whose n(r) peaks 
at a value of .1 ^m. The 8 a is the sum of all particles weighted by scattering effi- 
ciency peaks at .63 /xm, therefore, the 8 a is sensitive to a larger radius particle than 
CN measurements. Utilizing a Pearson correlation moment (Walpole and Meyers 
1985), the correlation for the two data sets is 0.56 for 59 samples. This positive 
correlation is a strong indicator of the validity of satellite retrieval techniques. Fur- 
ther examination of Fig. 4.8 and 4.9 reveal that the largest values for optical depth 
and condensation nuclei occur in the 10-50° N band, while the smallest values with 
the least variation occur south of 10° N. By separating the data into a North data 
set, the correlation coefficient was a positive correlation of 0.6 for 42 samples. For 
the South data set, there was a negative correlation of 0.44 for 17 samples. This 
negative correlation is indicative of the small data set with little variation and the 
differences between sensitivity of CN counts and 8 a measurements. 

The second largest feature for both optical depth and condensation nuclei is 
centered at 20° N. The major aerosol producing event in this region is the eruption 
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of the Kilauea volcano. Positive correlation between the two data sets verifies the 
satellite retrieval techniques. Contrary to the Northern region, the CN count is 
larger and the optical depth is smaller. Optical depth increases when the number or 
radius of the aerosol particles increase. The CN count increases when the number 
of aerosol particles increase. Therefore, it appears there are a greater number of 
small particles in the atmosphere. This may be a reflection of the settling of the 
low-level debris out of the atmosphere as it moves away from the volcano or it may 
be that the aerosol optical depth is made up of particles with larger radii at 170° W. 
Therefore, with a weaker source, more of the optical depth is a function of higher 
altitude particles. 

The peak at 10° N is most likely due to advection of continental sources by the 
prevailing NE trade winds (Bates and Charleson 1989). The change in optical depth 
and CN are three times smaller than in the westerly source off the Asian continent 
in the 40° N region. This appears to relate to the different type of continental source 
as well as to the effects of cloud processes in the tropical region. The peak near the 
Equator at 2° N is most likely biogenic in origin and is due to the upwelling, as there 
are no continental sources contributing in this region (Bates and Charlson 1989). 
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Figure 4.8: Optical Depth - Measured by NOAA-9 AVHRR within ± 3 

days of Ship’s Position 
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Figure 4.9: Condensation Nuclei - Measured by the R/V 
OCEANOGRAPHER 



34 



D. INTERCOMPARISON OF S 12 RATIO AND NSS-50]" 

Fig. 4.10 shows that the dominate features for the Su are maximums at 35° 
N, 19° N, 10° N, 1° N, and 2° S. Recall from earlier discussions that these elevated 
values indicate a greater ratio of small particles to large particles in the atmospheric 
distribution. Charlson, et al., (1987) describe a mechanism in which high oceanic 
DMS productivity by phytoplankton relates to increased numbers of small aerosol 
particles in the atmosphere. Therefore, one could infer that increased DMS pro- 
ductivity in a region (and the resultant increase in NSS-SC^ - ) should correlate 
with increased S 12 values. However, there are several factors which influence such 
a direct correlation. First, it has been shown that the primary sources for aerosols 
in the atmosphere for the North Pacific region (10-45° N) of this study are an- 
thropogenic, volcanic, or originate from atmospheric transport of Gobi desert dust 
storms. Therefore, biogenic activity in this region is either extremely limited or 
overshadowed by the continental sources. Second, the North Pacific Ocean region is 
noted for it’s reduced DMS productivity (Gross 1972). Third, during the cruise of 
the R/V OCEANOGRAPHER, it was noted that DMS production was at a winter 
low (Bates and Charlson 1989). Therefore, a direct correlation between DMS and 
S 12 is not seen in the North Pacific. 

In the equatorial region, the aerosol particles appear to be biogenic in origin. 
This is because the Inter-Tropical Convergence Zone (ITCZ) appears to be an effec- 
tive removal mechanism for continental sources of aerosols (Clarke 1989). Andrea 
(1985) found that DMS concentrations are very high in upwelled tropical waters. 
Additionally, Clarke (1987) observed that DMS production in the equatorial region 
is 3-5 times greater than in the North Pacific. Therefore, it is reasonable to expect 
that the NSS-^O^ - aerosols are biogenic in origin. 



35 



Fig. 4.11 shows the NSS-SO^ - data obtained from shipboard measurements 
on the R/V OCEANOGRAPHER. In the equatorial region, as expected, there is a 
maximum in sulfate measurements from 4° N to 2° S. The presence of this sulfate 
aerosol over such an extensive region remote from continental sources is consistent 
with production by biogenic sources such as DMS. The satellite derived Su has 
peak values at 1° N and 2° S and, in conjunction with the low optical depths for 
this area, from Fig. 4.8 provide strong support that these particles are biogenic in 
origin. It also coincides with Clarke (1987) in that there is an increase in atmospheric 
sulfur compounds in the equatorial region. Shipboard measurements of condensation 
nuclei also showed an increase in this region. This supports Charlson’s hypothesis 
concerning increased oceanic DMS production. 

The most prominent feature in Fig. 4.11 is the large peak in sulfate particles 
in the 20° N region. This maximum is a result of the volcanic emissions from the 
Kilauea volcano, and correlates with a similar maximum in 5n at 20° N. 

Shipboard sulfate measurements in the region of 40° N were limited due to the 
high winds and rough seas encountered. With the continental aerosol input from the 
westerly flow, it would be reasonable to expect higher NSS-S'O* - values. However, 
the aerosol particles in this region may not be sulfur based, which is indicative of 
atmospheric dust transport vice combustion sources from the Asian continent. Also, 
the 5 12 at 40° N is smaller, therefore, indicating that there are more large particles. 
Additionally, the maximum in S 12 located at 35 N may indicate that this region is 
anthropogenic in origin. The 38-44 N region maximums for both optical depth and 
condensation nuclei confirm the presence of a large number of aerosol particles in 
this region. The source of these particles is the westerly transport of continental and 
man-made aerosols due to the westerly synoptic flow pattern previously discussed. 
Interestingly, the peaks in optical depth at 42° N and 39° N probably represent 
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separate sources. The 42° N peak with a smaller optical depth than at 39° N has 
a larger APSI as shown in Fig. 4.10. The larger APSI with smaller optical depth 
indicates a greater number of smaller particles which would probably be of man- 
made origin. Conversely, the peak at 39° N has a larger optical depth with a smaller 
5 12 , therefore, the aerosols are larger and of a different source, most likely the Gobi 
desert “Kosa”. 

The last feature to discuss is the corresponding peaks at 10° N for sulfate 
measurements and .S 12 . The most likely source for these peaks is the NE tradewinds 
transporting continental sources of aerosols. The correlation of the above data points 
is another strong verification of the satellite retrieval methods. 
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Figure 4.10: Sn Ratio - Measured by NOAA-9 AVHRR within ± 3 days 

of Ship’s Position 
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Figure 4.11: NSS-SO^ Data - Measurements made by R/V 
OCEANOGRAPHER during April 1988 
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E. CHANNEL 1 AND CHANNEL 3 LOW CLOUD REFLECTANCE 



Fig. 4.12 is the channel 1 (0.63 /xm) low cloud reflectance data for the entire 
period. There is very little variation in the channel 1 reflectance values. Recall that 
the magnitude of the reflectance in this portion of the electromagnetic spectrum 
is heavily dependent upon cloud thickness and liquid water content (LWC) of the 
cloud, as well as the droplet size distribution. Therefore, these three variables must 
be offsetting each other, and little change in the reflectance is noted over the entire 
region. It is extremely difficult to obtain a correlation between reflectance and cloud 
droplet size in this channel. Fig. 4.13 and 4.14 are composite optical depths and 
S \2 measured by the satellite for the entire period. Higher optical depths and 5 12 
values are noted in the 10-40° N region. This would indicate that there are a greater 
number of aerosol particles with a smaller mean radius in the northern regions. Also, 
shipboard measurements of condensation nuclei increased in the northern region of 
this study. Channel 1 low cloud reflectance values do not show any effect of this 
variability. 

For channel 3 (3.7 /xm) reflectance, recall that the dependence on LWC and 
cloud thickness is not nearly as great as the dependence on droplet radius. Fig. 4.15 
is the low cloud reflectance of the entire region of this study. It is a significant result 
that low cloud reflectance is increasing from 10-40° N with a maximum reflectance 
value at 42° N. This is a strong indicator that changes in cloud reflectance based 
on droplet size distribution is more observable using the channel 3 information and 
strongly supports the satellite retrieval techniques developed by Mineart (1988). 

When CN (Fig. 4.9) is compared to the channel 3 low cloud reflectance, there 
is a direct correlation in the region of maximum values. Therefore, as the number of 
CN increase, the number of CCN increase, and result in higher channel 3 low cloud 
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reflectance. This is clearly demonstrated at 20° N, where the increase in CN due to 
volcanic activity results in channel 3 reflectivity increases. 
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Figure 4.12: Channel 1 Low Cloud Reflectance - Measured by NOAA-9 

AVHRR 
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Figure 4.13: Composite Optical Depth - Measured by NOAA-9 AVHRR 
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Figure 4.14: Composite S 12 - Measured by NOAA-9 AVHRR 
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Figure 4.15: Channel 3 Low Cloud Reflectance - Measured by NOAA-9 

AVHRR 
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V. CONCLUSIONS AND 
RECOMMENDATIONS 



This study utilized data from the RITS-88 cruise and NOAA-9 AVHRR to ver- 
ify satellite retrieval techniques, expand shipboard observations to a regional scale, 
and verify the relationship between aerosol characteristics and cloud reflectance 
properties. From these data, the aerosol characteristics were inferred and a better 
understanding of the radiative transfer properties of the atmosphere was obtained. 

Optical depth values show a positive correlation with the shipboard conden- 
sation nuclei counts. This provides strong justification for utilizing satellite derived 
optical depth measurements when estimating the radiative transfer characteristics 
of a region. When combined with the Si 2 values, these two techniques are very use- 
ful in determining the type and quantity of aerosols in the atmosphere. This would 
provide valuable information to global climate studies, as well as to the effects of 
aerosols on electro-optical systems. 

The relationship between solar reflectance and cloud droplet size distribution 
was illustrated. High channel 3 reflectivities are associated with high concentrations 
of aerosols, high concentrations of cloud droplets, and a shift toward smaller mean 
cloud droplet radius. This could provide a valuable tool for tracking the effects of 
major events, such as forest fires, industrial accidents, or volcanic emissions effect 
on the environment. 

The expansion of shipboard data to a regional scale via satellite measurements 
resulted in the identification of several strong aerosol sources in the North Pacific. In 
conjunction with this, it provided more evidence of the ability of the atmosphere to 
transport aerosol particles great distances. The volcanic plume from Kilauea volcano 
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was identified as a significant aerosol producing source for this region. The two 
major events of desert duststorms and volcanic emissions highlighted how a spatial 
view of the relative aerosol size distribution was possible from the S 12 and optical 
depth data. This made it relatively easy to distinguish between anthropogenic and 
naturally occurring aerosols. This would not have been possible utilizing only in- 
situ shipboard data and it highlights the valuable contribution of satellite retrieval 
methodology. 

Savoie and Prospero (1989) proposed that NSS-SO 4 - concentrations over the 
Pacific appear to be primarily controlled by biogenic sources in the ocean. Contrary 
to their conclusions, during the time period of this study, continental sources of 
aerosols had a greater impact in the North Pacific region. Continental influences 
would hinder any efforts to examine the relationship of DMS production to increasing 
the number of CCN in the atmosphere. Shema (1988) proposed that the equatorial 
region is a better place to study the relationships between DMS, aerosols, and clouds. 
Clarke (1989) proposed that the ITCZ appears to filter out continental sources of 
aerosols in the equatorial region. As a result, the aerosol particle population is 
mainly biogenic in origin, which are precursors to DMS production. Therefore, the 
study of Charlson’s biological feedback mechanism would be greatly enhanced in 
this region. 

The simple satellite retrieval techniques utilized in this study show that the 
effect of changing cloud reflectance due to spatial variations in aerosol concentration 
and composition can be incorporated into global climate studies using archived 
AVHRR data. 

Recommendations for future studies of this nature should compare in-situ 
measurements over the entire troposphere versus satellite data and incorporate CCN 
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measurements. This will give a better representation of the validity of satellite 
retrieval techniques. 
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